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^ ' Abstract 

1^ ■ Making use of the recently developed chiral power counting for the physics of nuclear matter [Tl[2] , 

ly^ , we evaluate the in-medium chiral quark condensate up to next-to-leading order for both symmetric 

^SJ ' nuclear matter and neutron matter. Our calculation includes the full in-medium iteration of the leading 

I — I , order local and one-pion exchange nucleon-nucleon interactions. Interestingly, we find a cancellation 

f^ ' between the contributions stemming from the quark mass dependence of the nucleon mass appearing in 

^^ , the in-medium nucleon-nucleon interactions. Only the contributions originating from the explicit quark 

mass dependence of the pion mass survive. This cancellation is the reason of previous observations 
concerning the dominant role of the long-range pion contributions and the suppression of short-range 
nucleon-nucleon interactions. We find that the linear density contribution to the in-medium chiral 
quark condensate is only slightly modified for pure neutron matter by the nucleon-nucleon interactions. 
C^ ' For symmetric nuclear matter the in-medium corrections are larger, although smaller compared to 

other approaches due to the full iteration of the lowest order nucleon-nucleon tree-level amplitudes. 
Our calculation satisfies the Hellmann-Feynman theorem to the order worked out. Also we address 
the problem of calculating the leading in- medium corrections to the pion decay constant. We find 
that there are no extra in-medium corrections that violate the Gell-Mann-Oakes-Renner relation up 
to next-to-leading order. 



1 Introduction 

The QCD ground state, or vacuum, is characterized by the presence of a strong condensate (0|gg|0) of 
scalar quark-antiquark pairs - the chiral quark condensate - which represents an order parameter for 
spontaneous chiral symmetry breaking in QCD. It is accepted on phenomenological grounds that the 
lightest pseudoscalar mesons, the pions, are identified with the Goldstone bosons of the spontaneously 
broken chiral symmetry [3]. Spontaneous chiral symmetry breaking results because the axial charge, 
Qa ~ J d^xAQ{x), does not annihilate the ground state. The coupling of the Goldstone bosons to 
the axial- vector charge is given in terms of the pion decay constant f-,^, which determines the chiral 
scale 47r/,r ~ 1 GeV that governs the size of the quantum corrections in Chiral Perturbation Theory 
(xPT) [4J. Chiral symmetry is also explicitly broken due to the small masses of the lightest quarks, u 
and d. An interesting issue for the discussion of the QCD phase diagram is the dependence of the chiral 
quark condensate on temperature and density. It is expected that with growing temperature, the quark 
condensate melts [5H3- There is an indication that this is also the case at zero temperature and increasing 
density as follows from the direct application of the Hellmann-Feynman theorem [8] to the energy density 
of a Fermi sea. A restoration of chiral symmetry is believed to be linked to a phase transition in QCD. 
A nonvanishing chiral quark condensate represents a sufficient condition for spontaneous breakdown of 
chiral symmetry. Note, that it is not a necessary one, such one would be the Goldstone boson (pion) 
decay constant, therefore chiral symmetry might as well be broken with a vanishing quark condensate. 
As long as the matrix element {Q\QWtt'^ (p)) = iSia{2TT)^ ftpoS{p) is not zero, the ground state is not 
left invariant by the action of the axial charge Q\ and spontaneous chiral symmetry breaking happens. 
In the previous equation |7r"(p)) denotes a pion state with Cartesian coordinate a, three-momentum p, 
energy po and ft is the temporal weak pion decay coupling. Many recent calculations in nuclear matter 
share the assumption that spontaneous chiral symmetry breaking still holds for finite density nuclear 
systems. The form of the chiral Lagrangians changes depending whether the chiral quark condensate 
{0,\uu + dd\Q) is large or small. In the former case we have standard xPT [91-lllj and in the latter 
generalized xPT J12pi3j should be employed. It has been shown that the first case holds [Hj in vacuum 
for SU(2) xPT. For modern applications of chiral symmetry to nuclear systems we refer to [15]. The 
in-medium chiral quark condensate for symmetric nuclear matter in the linear approximation |16fl21j 
decreases with density as 1 — (0.35 it 0.09)p/p0) with the error governed by that of the knowledge in the 
pion-nucleon sigma term, a = mdm/drh = 45 it 8 MeV |22y23j. with m = (rn-^ + m(])/2 the mean of the 
u and d quark masses, m the nucleon mass and po — 0-16 fm~^ the nuclear matter saturation density. 
Thus, if the low energy linear decrease in the quark condensate is extrapolated to higher densities the 
quark condensate would vanish for p = (2.9 it 0.7)po and standard xPT would not be appropriate. Of 
course, higher order corrections could spoil this tendency. Hence, the calculation of the in-medium quark 
condensate is of great interest beyond the linear approximation. The quark condensate in the nuclear 
medium and the quest for restoration of chiral symmetry at finite baryon density (and temperature) has 
a long and outstanding history, see e.g. [T6H2T|[24H32] . 

In this work we concentrate first on the application of in-medium baryon xPT to the calculation of the 
quark condensate and go beyond the linear density approximation. However, our approach [Ij compared 
with previous works offers two novel features: i) It follows a strict chiral power counting that takes into 
account both long- and short-range (multi-)nucleon interactions. It is applicable both in vacuum and 
in the medium so that a clear connection between the two cases is established and used, ii) We do not 
take as starting point the Hellmann-Feynman theorem but directly apply the power counting mentioned 
to the problem of the calculation of the quark condensate. The fulfilment of the Hellmann-Feynman 
theorem is a consequence of the consistency of the approach order by order. In this way, the dependence 
on the quark mass of the input parameters in the theory, like /jr, Qa, nucleon and pion masses, etc, is 



built in. This power counting has been previously applied to the problems of the pion self-energy in 
the nuclear medium and the equation of state of nuclear matter. For the former issue [1] we found a 
cancellation of the next-to-leading order (NLO) contributions, including those stemming from the in- 
medium nucleon-nucleon interactions. This suppression is interesting since it allows to understand from 
first principles the phenomenological success of fitting data on pionic atoms with only meson-baryon 
interactions [33H35] . Regarding the equation of state of nuclear matter it is interesting to point out that 
saturation for symmetric nuclear matter and repulsion for neutron matter is obtained in ref. [2] at NLO. 
The resulting curves of the energy per particle as a function of density for neutron matter agrees with 
sophisticated many-body calculations ^36j in terms of a subtraction constant go that acquires its natural 
expected value, close to —0.5 m^. For the case of symmetric nuclear matter the experimental properties 
of the minimum for the bounding energy are reproduced fitting go ~ —m-^, while the analogous parameter 
go keeps its just mentioned natural value. As a result of the present and the previous applications p] f2ll21] 
the calculation of the in-medium corrections to the coupling of the pion to the axial- vector current up to 
NLO is straightforward. As we show below, the linear reduction with density of ft as obtained in ref. |21] 
holds up to NLO. 

After this introduction, we present the power counting formula used and the Feynman diagrams 
contributing to the in- medium chiral quark condensate to NLO in section [2j The contributions stemming 
from pion-nucleon chiral dynamics are the object of section [3l The terms due to the nucleon-nucleon 
interactions are calculated in section [H Results and discussions thereof are given in section [5l The in- 
medium pion decay constant and the Gell-Mann-Oakes-Renner relation are discussed in section [6l The 
last section is dedicated to offer our conclusions. 

2 In-medium chiral power counting and diagrams 

For determining the set of diagrams to be calculated in the evaluation of the quark condensate in the nu- 
clear medium up-to-and-including NLO contributions we make use of the chiral power counting developed 
in [5 

K V 

v = A-E + ^{ii + ni - 4) + ^{di + Vi + uJi-2) + Vp . (2.1) 

i=l i=l 

Li this way, E is the total number of external pion lines and V^, V and Vp are the meson- meson, meson- 
baryon and in-medium generalized vertices |37j . in this order. Li simple terms, an in- medium generalized 
vertex corresponds to a closed nucleon loop that could contain an arbitrary number of meson-baryon 
vertices inserted. In the first summation on the right-hand-side of the previous equation the symbols ii 
and Hi are the chiral order and number of pionic lines of the i meson-meson vertex, respectively. In the 
sum over the meson-baryon vertices, di is the chiral dimension of the i^^ meson-baryon vertex, Vi is the 
total number of mesons lines attached to it, including both pion and heavy meson lines, while Wj is the 
number of the latter ones only. The heavy meson lines correspond to auxiliary fields responsible for the 
local multi- nucleon interactions when taking their masses to infinity [Tl[38]. 

Eq. ()2.ip counts every nucleon propagator as 0{p~'^) instead of 0{p~^), as in the standard counting 
used in Baryon xPT |391I40|. In this way, the infrared enhancements associated to the large nucleon 
mass are taken into account from the onset. Despite baryon propagators are counted as 0{p~'^) it is 
important to stress that eq. (|2.ip is bounded from below [2]. According to eq. ()2.ip the number of lines 
in a diagram can be augmented without increasing the chiral power [l] by adding i) pionic lines attached 
to lowest order mesonic vertices, ii = rii = 2, ii) pionic lines attached to lowest order meson-baryon 
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Figure 1: Contributions to the in-medium chiral quark condensate up to NLO or 0{p^). The scalar source with 
zero momentum is indicated by the wavy hne and pions by the dashed ones. A wiggly hne corresponds to the 
nucleon-nucleon interaction kernel (given in fig. |3]) whose iteration is denoted by the ellipsis. 



vertices, di 



1 and iii) heavy mesonic lines attached to lowest order bihnear vertices, di = 0, 



Wj = 1. In this way, resummations of infinite strings of diagrams are required, and these resummations 
signal the appearance of non-perturbative physics. As a major difference between this counting and the 
Weinberg one |41ll42j . the former applies directly to the physical amplitudes while the latter does only to 
the potential. It is also important to stress that z/ in eq. (j2.ip is a lower bound for the actual chiral power 
of a diagram, ^u, with ^ > y. The chiral order might be higher than v because nucleon propagators, 
always counted as 0{p^'^) in eq. (j2.ip . could actually count as ©(p^^) in some cases. If this is so, every 
0{p~^) baryon propagator increases the actual chiral order by one unit as compared with v eq. (12. ip . 
The important point is that all the diagrams with a chiral order ^ are a subset of the diagrams with 
u < ^. We refer to [l] for more details on eq. (j2.ip . including its derivation. 

The explicit chiral symmetry breaking due to the quark masses is incorporated by s{x) = A4. + 5s[x) 
with M = diag(mu,md) the quark mass matrix. From the generating functional in the presence of 
external sources, Z{v,a,s,p), the quark condensate is obtained by partial functional differentiation 



{n\qiqj\n) 



5sij {x) 



Z{v,a,s,p%, 



a,s,p=0 



(2.2) 



with qi = u and q2 = d quarks. Notice that the quark condensate has a global minus sign compared 
to diagrams calculated by using ordinary Feynman rules. See ref. |37] for a derivation of the generating 
functional Z{v,a,s,p) in the nuclear medium making use of functional methods, although keeping only 
pion-nucleon interactions. We also refer to [21] for more details on the use of external sources in relation 
with in- medium xPT calculations. 



The quantity we will calculate in the following is the in-medium correction H to the chiral quark 
condensate given by 

mqmqiqj\n) = mg{0\qiqj\0) + niqE , (2.3) 

where niq is the mass of a certain quark flavor. Furthermore, the nuclear matter ground state is denoted 
by |r2), to distinguish it from the vacuum ground state |0). The resulting contributions are denoted by 
Hj and are shown in fig. [TJ To determine the set of diagrams needed for the calculation of the in-medium 
chiral quark condensate up to NLO we proceed by increasing Vp step by step in eq. (12. ip . For each Vp 
we then determine the possible configurations of vertices and lines according to eq. (j2.ip . The resulting 
diagrams are shown in fig. [TJ They can easily be identified by considering the corresponding vacuum 
diagrams. For Vp = 1 the first diagram is the lowest-order contribution to the nucleon sigma-term. The 
rest of diagrams in fig. [1] are NLO. The diagrams 2 and 3 stem from the one-pion loop self-energy with 
the scalar source attached to a nucleon or pion propagator, respectively. For Vp = 2 one can conveniently 
think of the nucleon-nucleon scattering in the presence of a scalar source in vacuum. Then proceed by 
closing the diagrams, which corresponds to sum over the states in the Fermi seas of the nucleons. The 
upper row of diagrams 4, 5 and 6 in fig. [1] involves the direct nucleon-nucleon interactions while the 
lower row originates from the exchange part. In the following we indicate by Hj the contribution to 
the in-medium quark condensate due to the diagram i in fig. [TJ One could think of similar Feynman 
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Figure 2: These diagrams are zero because of the antisymmetric structure of the Weinberg- Tomozawa vertex 
coupling in the pion indices. For more details see the text. 

graphs to diagrams 2 and 3 of fig. [TJ but involving two-pions in one vertex from the Weinberg- Tomozawa 
term of C^p^ eq. (12.61) . These are depicted in fig. [2j However, these diagrams are zero because of the 
antisymmetric isospin structure of the Weinberg- Tomozawa vertex, proportional to £imcT^- The same 
type of pion is involved in the tadpole loop. This is clear for the diagram a) but also for diagram b) due 
to the diagonal structure of the vertex coupling the scalar source with two pions, eq. (|2.10p . As a result, 
both diagrams are zero when the pion indices are contracted with the antisymmetric tensor. 
The in-medium nucleon propagator reads 



with Ei\i) = k^/2m, the non-relativistic nucleon kinetic energy. In this equation the subscript i^, is the 
third component of the isospin of the nucleon, so that i^, = -1-1/2 corresponds to the proton and is = —1/2 
to the neutron. The symbol ^jg is the Fermi momentum of the Fermi sea for the pertinent nucleon. We 
consider that in vacuum isospin symmetry is conserved, so that all the (vacuum) nucleon and pion masses 
are equal. The proton and neutron propagators can be combined in a common expression 

Go(fc) = j;(^ + ^3T3)Go(fe).3. (2.5) 



We make use of the Heavy Baryon xPT (HBxPT) pion-nucleon Lagrangian of |44j up to 0{p 



C^^^ - N 



'-'■kN — ^^ 



iDq - —a ■ u 



N (2.6) 



where (. . .) denotes the trace in isospin space and the ehipses represent terms that are not needed here. 
The covariant derivative is given by D^ = S^+F^, where the chiral connection is V^ = 2 [m^, 9/itt] — §it^(i'^+ 
o-^)u— ^u{Vfj, — a^)v) . The chiral vielbein is given by u^ = i{v) ,d^u} + v) [v^ + a^)u — u{v^ — a^)v) . The 
couphng of a scalar source to pions is also required. At lowest order in the chiral counting this coupling 
arises from 



6^) = fl{^^) + ... (2.7) 

The matrix field of scalar sources s{x) = ||sij(x)||, is included in the operator x+ [10] - 

X+ = u^XU^ + ux^u 
X = 'iBsix) . (2.8) 

The parameter B is related to the strength of the vacuum chiral quark condensate and the weak pion 
decay constant fj^, both in the chiral limit, via 

B6,j = -{0\mj\0)/f^. (2.9) 

The vertex coupling the scalar source Sij to two pions of Cartesian coordinates a and b can be evaluated 

(2) 

straightforwardly from —ll\i^ eq. ()2.7p . with the result 

Bblrr^'' ■ (2.10) 

This vertex is diagonal both in the scalar source as well as in the pion fields. 

The nucleon-nucleon scattering amplitude at lowest order, ©(p*^), is needed to match with our aim 
of calculating the NLO contributions to the chiral quark condensate in the nuclear medium. Note that 
the nucleon-nucleon scattering amplitudes first appear in diagrams with V^ = 2 so that it starts to give 
contributions to the quark condensate already at NLO. These amplitudes originate from the Lagrangian 
with four nucleons, without quark masses or derivatives |42) 

^SJv = -\cs{NN)[NN) - hjT{NaN)(NaN) , (2.11) 

and from the one-pion exchange with the lowest order pion-nucleon coupling. Its sum is represented 
diagrammatically in the following by the exchange of a wiggly line as in fig. [3l For explicit expressions 
of these amplitudes see [2] , where we also compare with vacuum nucleon-nucleon scattering data. 

3 Vp = 1 contributions 

First we consider the contributions to the in-medium chiral quark condensate from pion-nucleon chiral 
dynamics. They are depicted in diagrams 1-3 of fig. [TJ For the evaluation of the different diagrams we 
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Figure 3: Leading order interaction kernel: the exchange of a wiggly line between two nucleons indicates in the 
following the sum of the 0{p'^) local and the one-pion exchange contributions. 

need the vertex with the scalar source Sij coupling to a pair of nucleons I, m. It can be readily worked 

(2) 

out from —C}^lf eq. (|2.6p . with the result 

-2B[2ci5ij5lm + C^Tji ■ Tim\ , (3.1) 

where we have taken into account that 25ii5jm — ^ij^im. = ^ji ■ ^im- The diagram 1 then yields 

Hi = -2B [2cA,ipp + Pn) + C5(T3)y(/9p " Pn)] = ^1 + Sf , (3.2) 

where the pp and /?„ are the proton and neutron densities given by Pp{n) = Cp{-„)/37r^. Notice that the 
isospin breaking contribution proportional to C5 only involves the Pauli matrix rs, so that for i ^ j the 
contribution vanishes, as required. The contribution in eq. ()3.2p proportional to ci is denoted in the 
following by r}-^ and that proportional to C5 by r.f . The superscripts is and iv refer to the isoscalar and 
isovector character of these contributions, respectively. 

Let us proceed to the evaluation of the NLO diagrams 2 and 3 of fig. [TJ Diagram 2 originates by 
dressing the in- medium nucleon propagator with the one-pion loop nucleon self-energy, S'^. It is given 

by 

H2 = -2iB I -0^6'^'^ Tr { [2ci6ij + c^Tji ■ r] G,{k) S^ Go{k)] = Ef + E^^ , (3.3) 

with the convergence factor e**^ '', r/ — )• 0+, associated with any closed loop made up by a single nucleon 
line |43) . The trace, indicated by Tr, acts both in spin and isospin spaces. Similarly as in eq. (13. 2p the 
term proportional to ci is denoted by Hg* and that proportional to C5 by H™. The contribution from 
diagram 2 with only free parts in all the nucleon propagators involved, including that in Ti^ , vanishes. 
This is discussed in detail in ref. [2] for similar diagrams that appear in the calculation of the in-medium 
pion self-energy and nuclear matter energy. Briefly, it follows just by closing the integration contour 
on the complex k^ half-plane opposite to that where the poles lie. Another important point to keep in 
mind is that the contributions with only Fermi sea insertions in all nucleon propagators is part of the 
Vp = 2 contribution of the crossed exchange part of diagram 4 in fig. [TJ This is shown diagrammatically 
in fig. 3 of [2j, where the two external pion sources should be replaced by the scalar source for the case at 
hand (similarly, see fig. [H for diagram 4.) The different Vp = 2 contributions are evaluated in section HI 
Consequently, only those terms that involve simultaneously free-space as well as density-dependent parts 
of the nucleon propagators are considered here. For the isovector part the calculation of the contribution 
lis in ref. [2\ applies straightforwardly, just by replacing the vertex function. The same procedure applied 
to the isoscalar part drives to its cancellation, E'if = 0. For the same reason as in |[2| Sg^ turns to be of 
0{p^) or N'^LO. This is due to the appearance of the derivative of the free part of the nucleon one-pion 
loop, ST, with respect to energy. This derivative is finally suppressed by one chiral order [2]. 

For diagram 3 in fig. [1] the scalar source is attached to the pion propagator involved in the one-pion 
loop nucleon self-energy. In the same way as discussed above for diagram 2, the contribution with only 
the free part for all nucleon propagators vanishes, while that with all nucleon propagators as in-medium 
insertions are accounted for by the exchange part of diagram 6 of fig. [H see fig. [H This is calculated 






Figure 4: The equivalence between diagram 3 and the crossed part of the one-pion exchange reduction of diagram 6 
of fig. [His shown. The diagram in the middle is an intermediate step in the continuous transformation of the diagram 
on the left-hand-side to the one on the right-hand-side. 



in section [H Then, we consider here the part of diagram 3 having one free nucleon propagator and the 
density dependent part of the other. Let us recah that the couphng of a scalar source to two pions was 
already given in eq. (j2.10p . There is a close relationship between diagram 3 of fig. [1] and the contribution to 
the nuclear matter energy from the one-pion loop nucleon self-energy, diagram 2 in fig. 17 of [2], denoted 
by E^- Of course, this is a requirement from the Hellmann-Feynman theorem. It is straightforward to 
check that 



Ha = iB6i 



u 



d k_^if,o^^^ 
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Goik) 



dm'i 



B5i 



d£2 



(3.4) 



with the derivative affecting only the explicit dependence of ST on the pion propagator, and not including 
the implicit one from the nucleon mass dependence on it. One has to subtract the value of the one-pion 
loop nucleon self-energy at fc'' = since we are using the physical nucleon mass. After performing the k^ 
integration, one has the expression 



where to = -©(k) and the partial derivative is given by 
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(3.6) 



with h = m^ — bj'^ — ie and e — )• 0+. Notice that H3 is an 0{p^) or N^LO contribution because dTi^Jdrn^ = 

Oip^). The 0{j>) contribution from the term —?>TTm.,^ is cancelled by that coming from 37rv6. We then 
conclude that the only contribution with Vp = 1 up-to-and-including NLO is given by Hi, eq. (13. 2p . 



4 Vp = 2 contributions 

We now consider those next-to-leading order contributions to the in-medium chiral quark condensate 
that involve the nucleon- nucleon interactions. They are depicted in diagrams 4-6 of the last two rows of 

fig.m 



4.1 Contributions S4 and S5 

Diagrams 4 and 5 are analogous to diagrams a) and c) of fig. 4 in [i\. There it has been shown that those 
diagrams cancel each other. This argument was developed in [2] within a specific method for resumming 
the iteration of the wiggly lines for the non-perturbative nucleon-nucleon interactions. We expect that 
this cancellation also takes place for the case of the in-medium chiral quark condensate. However, the 



vertex coupling two nucleons with the scalar source, eq. ()3.ip . has both an isoscalar and an isovector 
term while for the aforementioned cancellation in the case of the pion self-energy only an isovector vertex 
was involved. Here we want to show on general grounds that the cancellation also takes place for the 
in- medium chiral quark condensate, following similar arguments as those in [T]. 
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b) 



Figure 5: Contribution to the chiral quark condensate with a two-nucleon reducible loop. The scalar source 
couples outside the loop for diagram a) and inside it for diagram b) . 

Diagram 4 of fig. [1] can be written in terms of the nucleon self-energy due to the in-medium nucleon- 
nucleon scattering, denoted by S . It reads 



-2iB / - — —e 

(27r)4 
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^ Tt{{2ci6i, + C5Tji ■ f)Go(A;i)S^^Go(fci)} 



The expression for T,ai,NN, corresponding to the self-energy of a nucleon with isospin ai, is [T] 



^ai,NN 
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02,0"2 
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(4.1) 



(4.2) 



Here Q2, (T2 correspond to the running third components of isospin and spin, respectively, while ai is 
the third component of the spin of the external nucleon. In addition, T^lali^'^^^^) refers to the elastic 
in-medium nucleon-nucleon scattering amplitude for Nai,cTiiki)Na2,a2{k2) — ^ -^ai,(Ti(^i)-^a2,o-2(^2)- We 
also make use of the identity dGo{k)/dk^ = —Go{k)'^. As a result, eq. (|4.1|) can be expressed as 
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(4.3) 



where in the last equation an integration by parts in k^ has been performed. In order to see the cancella- 
tion between diagrams 4 and 5 of fig. [Diet us proceed similarly as in ref. [Ij, for the case of the in-medium 
pion self-energy, and consider first the case with only one reducible two-nucleon diagram, fig. [5l The 
contribution of diagram a), H4, can be readily worked from eq. (j4.3p by reducing T^^^l{ki,k2) to its 
one-loop calculation. It results 



^4 



2^ E E /^^e^'°''e^'°''G'o(A:i),,(2ci5.,+C5(r3).,(r3)„,„jGo(fc2).2 



01,02 o"l,a'2 



d 



E 



_dS_ 



K.ia2;aio^(9)G'o(A;i - q)a[Go{k2 + ?)a^ V^/^a^jai 



02 > 



(4.4) 



where Vap-^^yS corresponds to the wiggly hne with the first pair of labels belonging to the outgoing nucleons 
and the second pair to the in-going ones. Note that a symmetry factor 1/2 is included because V contains 
both the direct and exchange terms. In addition, V also will depend generally on spin. Both isospin and 
spin indices are globally indicated in the labels of V with Greek letters. The quantity in squared brackets 
in the previous equation is T^lal^^ii^^) at the one-loop level. It is not necessary to consider T^l^l at 
the tree-level, consisting of the exchange of one wiggly line, fig. [3l because it is then independent on fcj* 
so that the derivative with respect to k^ is zero. 

One can similarly write down the contribution from the diagram b) of fig. [5l denoted by Sg . It reads 
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X {2ci5ij + C5{T3)ij{Ts)a'^a[)Go{ki - q)a[Go{k2 + q)a'^Va'^ 



a'„;aia2 ^ 



(4.5) 



with the global symmetry factor 1/2 from closing the lines. The appearance of the derivatives with respect 
to ki is again due to the fact that the propagator attached to the scalar source appears squared. The 
loop integrals between the squared brackets in eqs. (j4.4p and (|4.5p are typically divergent. Nevertheless, 
the parametric derivative with respect to k^ can be extracted out of the integral of eq. (j4.5p as soon as 
it is regularized. Summing over all isospin states makes clear that the position of the vertex associated 
with the coupling of the scalar source to two nucleons, either inside or outside the squared brackets, 
does not yield a difference. For that one has to keep in mind that T^^^'^{ki,k2) = T^^^^{k2,ki), due to 
the Fermi-Dirac statistics for a pair of two-nucleon states. Furthermore, notice that all the indices and 
four-momenta associated with the nucleons 1 and 2 are summed and integrated, respectively. In this way 
eqs. ()4.4p and ()4.5p mutually cancel|*l 






Figure 6: After performing the integration by parts in eq. (|4.3p the derivative with respect to k^ acts onto the 
scattering amplitude. This gives a sum of derivatives acting on two-nucleon reducible loops, indicated by the 
crosses. When the derivative acts on a baryon propagator the latter becomes squared. In this way, the first 
diagram in the second row of the figure equals the one of the first row but with opposite sign and they cancel each 
other. The same applies to the second diagrams in both rows. 

This process of mutual cancellation can be generalized to any number of two-nucleon reducible loops, 
using the same argument as given in ^ for the case of the in-medium pion self-energy. An n + 1 iterated 



*^ Alternatively, note that the term proportional to ci in eq. (|4.5[l is the same independently of its location because is 
proportional to the identity matrix. That part proportional to C5 is purely isovector, as those terms stemming from the 
Weinberg- Tomozawa vertex in refs. [Ull]- As shown there only the states with I3 = ±1 give contribution so that C£,{t3)^i ^1 
in eq. (|4.5I) gives rise to C5/3, which then trivially factorizes for the separate contributions of the proton-proton and neutron- 
neutron states. 
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wiggly line exchange implies n two-nucleon reducible loops. The scalar source can be attached to any 
of them for H5, while for H4 the derivative with respect to k^ can also act on any of the loops. This is 
exemplified in fig. [6] for the case with two two-nucleon reducible loops. Hence, 

H4 + H5 = . (4.6) 

The basic simple reason for such cancellation is that while for H5 there is a derivative acting onto the 
scattering amplitude, H4 involves an integration by parts in order to do so, which then introduces an 
extra minus sign. 

The previous cancellation is also explicitly obtained making use of Unitary xPT applied to nuclear 
matter as developed in ref. \2\. The intermediate result of [2j can be used straightforwardly by considering 
that instead of the Weinberg- Tomozawa vertex and pion-nucleon Born terms, used in the problem of the 
in-medium pion self-energy, one has the nucleon vertex of eq. (j3.ip with the scalar source. The partial 
wave decomposition of H4 reads 

J,i,S,Iai,a2 -^ *- ^'' ^ ^^ 

X GoihUGo{k2)a,m -^^^^ ' ^ , (4.7) 

with ai + a2 = -^3, the third component of the total isospin / of the two-nucleon state. Other symbols 
used are J the total angular momentum, i the orbital angular momentum and S the total spin of the 
nucleon-nucleon pair. We also employ the kinematic variable A = 2ma^ — a? , where a = {ki-\-k2)/2 and a 
is the three- vector made by the spatial components of a. For on-shell scattering A = p^, with p the three- 
momentum in the two-nucleon rest frame. For the term proportional to C5 due to its isovector character 
only the difference between the proton-proton and neutron-neutron contributions survives (Is = for 
proton-neutron pairs). For the expressions at 0{p^) and Oi^p"^) we plug into eq. (|4.7p the derivative 
dTjj/dA at leading and next-to-leading order, according to eqs. (5.24) and (5.25) of [2], respectively. For 
the contribution H5 the general partial wave decomposition yields 

"^ Ak-!r) AkJ.n 1^ ri^ \ /~i i -i^ 

^2)a2 

J,l,SJai,a2 -^ ^^^^ '^^^'' 

(4.8) 



-5 = -^ E E (2^ + 1)X(^^^)^ / (tH(lt4 ^^'°^^^'^^go(fci)..Go(fe2 












The expressions for Dfj at LO and NLO are given in eq. (5.25) of [2] and ^fj is calculated order by order. 
For that one has to work out the partial wave decomposition of the two-nucleon reducible diagram with 
the scalar source attached to one of the nucleon propagators inside the loop. Following those calculations 
one obtains for the leading order (LO) and next-to-leading order (NLO) 
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mdLiQ 



DLw = -4B[2cAj+hc5iT3) 



«jj 



dA 



Dlfj = -4B [2ci5,j + /3C5(r3).j] ^^^^ . (4.9) 

The contribution H5 at 0{p^) and 0{p'^) results when the just given expressions for ^j/, in order, are 
inserted in eq. ()4.8p . It is straightforward to check that they exactly cancel with H4 accordingly. 
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4.2 Contribution Sg 

Let us now consider the calculation of diagram 6 of fig. [H where the scalar source is attached to an 
exchanged wiggly line. Since the scalar source coupling to a local term is of higher order, it only couples 
to the pion exchange lines here. The appearance of the pion propagator squared with a zero momentum 
scalar source leads to its derivative with respect to the pion mass. The local term is independent of 
the pion mass, so we may act the derivative onto the explicit dependence on the pion-mass of the full 
nucleon-nucleon scattering amplitude. The implicit dependence of the latter on m^ due to that of the 
nucleon mass is the content of the diagrams 4 and 5 of fig. [H which have been shown above to cancel 
mutually. We can write 

=6 - ^B6^, 1^ 2^ 7^7^" ^ ^''Go(/ci)„,Go(fc2)a. Q^^ B6,,^^ . (4.10) 

rvi rvo (T-i /TO .^ V / V / TT TT 



Ql,Q2 0'1,0"2 



Where £^ is the contribution to the nuclear matter energy due to the in-medium nucleon-nucleon inter- 
action at leading order, evaluated in ref. p]. In order to proceed we have to evaluate the derivative of 
the nucleon-nucleon scattering amplitude with respect to m^. We decompose the in- medium nucleon- 
nucleon scattering amplitude in a sum over nucleon-nucleon partial waves Tji{£',i, S), as previously done 
for H4, eq. ()4.7p . We distinguish now between i and i' , corresponding to the the initial and final orbital 
angular momentum, in this order. In terms of the nucleon-nucleon interaction kernel Njj{i',i,S), the 
nucleon-nucleon partial wave Tjj in Unitary xPT is given by [2] 

Tjj{£',e,S) = [1 + Nji{e',£,S) ■ Lio]'^ ■ Nji{£',i,S) = Dj} ■ Nji , (4.11) 

with Lio the nucleon-nucleon unitarity scalar function |2J. This function does not depend explicitly on 
the pion mass. We have also introduced Djj = 1 -|- Njj ■ Liq. Eq. (j4.1ip can also be rewritten as 

Tji = Nji - Nji ■ Lio • Tjj . (4.12) 

Taking the derivative with respect to the explicit dependence on m^ on both sides of the previous 
expression leads to 

f4 = i^^ir' ■ f^ • [Dji]-' ■ (4.13) 

In terms of this equation the partial wave decomposition of Hg at NLO reads 



(fh d^k^ ii,o„ ii.o„„ ,. , „ ,, . dTjj 



J,e,S,Iai,a2 J \ ) y ) TT 



(4.14) 



LO 



On the problem of evaluating and regularizing this amplitude we refer to the related detailed discussion 
in [2] for the calculation of £3. 

5 Discussion and results 

It follows from our calculation for the in-medium corrections shown in fig. [1] that the chiral quark con- 
densate up-to-and-including NLO is given by 

mg{Q\qiqj\Q) = mq{0\qiqj\0) + mg{Ei + Eq) . (5.1) 
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Here we have taken into account that H2 and H3 are indeed one order higher (NNLO), while H4 + H5 = 0. 
Both Hi and Hg are clearly connected with the corresponding contribution to the nuclear matter energy, 
as required by the Hellmann-Feynman theorem 

mq{n\qiqj\n) - mg{0\qiqj\0) = 



?(*«i + (^»)«i)('""+^'' f^-^) 



with m = {niu — m^j^ and E the energy density of the nuclear matter system. This constraint is fulfilled 



in our case, where Hi, eq. (|3.2p . is the leading derivative with respect to m^ of the nucleon mass^J In 
turn, Hg corresponds to the explicit derivative of the nuclear matter energy density due to the nucleon- 
nucleon interactions with respect to m^, eq. ()4.10p . Notice, that the implicit dependence on m^ of the 
in-medium nucleon-nucleon interactions does not give any contribution to H because of the discussed 
above mutual cancellation between H4 and H5 at NLO. 

The term 'E}^ can be written directly in terms of the pion-nucleon a term with the result 



0"(Pp + Pn 
= 1 = -\U|C/i(/j|U; — -2-2" 

i'K"''T! 



(0|g.g.|0) T2 2 ' (5-3) 



using the Gell-Mann-Oakes-Renner relation (GMOR) [IS], rri^f^ = — 2m(0|gigj|0), valid at lowest order 
in the chiral expansion, see e.g. |47] . Regarding the numerical value of a one has the earlier extraction 
a = 45 lb 8 MeV of |221l23j or the one from the more recent partial wave analysis of pion-nucleon 
scattering [18] cr = 64 ± 7 MeV. In both cases the dispersion analysis of the nucleon scalar form factor 
of |22j to estimate the departure between the sigma-term and the vrA^ scattering amplitude evaluated 
at the (unphysical) Cheng-Dashen point is used. In asymmetric nuclear matter there are additional 
isospin symmetry breaking terms proportional to C5 stemming from diagrams 1 and 2, giving rise to 
the NLO contribution H™, eq. (j3.2p . This contribution, which distinguishes between the uu and dd 
quark condensates, are suppressed because they are proportional to the numerically small low-energy 
constant C5 = —0.09 it 0.01 GeV~^ [39]. For symmetric nuclear matter this contribution is zero since it 
is proportional to the difference of proton and neutron densities. 

The energy density and energy per particle, £/p, have been calculated in pj. In this reference, 
saturation of symmetric nuclear matter is obtained from the consistent chiral power counting, eq. (j2.ip . 
In addition, a remarkably good agreement with sophisticated many-body calculations [50] was obtained 
for the equation of state for neutron matter in terms of just one free parameter, a subtraction constant 
called go- The appearance of this subtraction constant is due to the purely free part of the scalar two- 
nucleon reducible loop, Lio already introduced in eq. (j4.1ip . Its free part, denoted by g{A), is a divergent 
function which requires a subtraction [2j. Taking D as the subtraction point one has 

g{A) = g{D) - \, / dk" 



4^2 J^ (fc2 _ ^ _ ie) (fc2 _ ^ _ ^g) 

,^, im / r- . /— — \ ,m\/A . ,. 

= aiP)--[^^-^^\)=g,-^^^. (5.4) 

(fa, eq. (7.36) in [2], depends on g^ both implicitly, due to the dependence of the partial wave amplitudes 
Tji{i',i, S) on go, cf.eq. ()4.1ip . and explicitly, in a linear manner. For the latter dependence the symbol 
5o was introduced in ref. [2j (to which we refer for a more detailed discussion on gQ and ^o)- For pure 
neutron matter both go and ^0 are obtained with the value go = go — —0.6 m^ , very close to their expected 



* At lowest order in the chiral expansion m^ = 2Brh so that in eq. (|5.2p one can make use of the operator m^d/dm^ 
instead of rhd/drh. 
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Figure 7: (Color online.) £ / p for neutron matter (left panel) and for symmetric nuclear matter (right panel). The 
two (magenta) dotted lines correspond to ref. |3^ and the solid ones are our calculation. 

natural size qq ~ —0.5 'm~'^. The case of symmetric nuclear matter requires some fine-tuning of qq, with 
a final value of qq ~ — m^, while Tjq is kept on its expected value of —0.5 m^. With that, results in perfect 
agreement with experiments were obtained for the saturation density, energy per particle and compression 
modulus p] , see fig. [7| where the equations of state for neutron (left panel) and symmetric nuclear matter 
(right panel) are shown. The magenta lines are from the many-body calculations of refs. [36p50j employing 
to so-called realistic nucleon-nucleon potentials. The agreement between these references and our results 
is remarkable. We use our previous evaluation of £^3 for calculating Hg, eq. ()4.10p . by taking the derivative 
with respect to the explicit dependence of E^ on m^. It is important to stress that the in-medium quark 
condensate does not depend on qq because this parameter multiplies a quantity which is independent of 
the pion mass, denoted by Sqo^ in ref. [2]. The latter originates from the limit p^ — )• 00 of the one-pion- 
exchange tree level partial waves at lowest order and this limit is independent of m^. Making use again 
of the GMOR relation, Hg eq. ()4.10p can be written as 



1 dSj 



=6 = -{0\qiqj\0) 



Putting together the previous equation with Ef, eq. (j5.3 
in-medium and vacuum quark condensate to NLO reads 



a + Oiml)) 



(5.5) 



and H™, eq. (13. 2|) . the ratio between the 



{^\qiqj\n) 
(Ol^i^ilO) 



1 



Cr{Pp + Pn) , '^C^iTs)ij{pp- Pn) 1 dSs 



fnl^l 



+ 






f^dml 



(5.6) 



The input parameters are qa = 1-26, /^ = 92.4 MeV, m^ = 138 MeV, tun = 939 MeV and a = 
45 MeV. In addition the values of the subtraction constants go that reproduce better the equation of 
state of symmetric nuclear matter and neutron matter are used. These values are [2] (70 = —0.97 and 
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Figure 8: (Color online.) The ratio between the in-medium and vacuum chiral quark condensate, {Q\qq\0) / {Q\qq\Q) , 
eq. (j5.6p . for neutron matter, left panel, and symmetric nuclear matter, right panel. Left panel: The (red) solid and 
(cyan) dot-dashed lines are our ftdl results for (l^j-uuiri) and {n\dd\il), respectively. The (black) dashed and (blue) 
dotted lines correspond in the same order to the linear density approximation (Si) only. Right Panel: The (red) 
solid and (black) dashed lines are the full results (with go = —0.97 m^) and the linear approximation, respectively. 
The (blue) dot-dashed line is the calculation for go = —0.5 m^. All the curves shown employ a = 45 MeV [22 . 



go = —0.62 m~^, respectively. In fig.[8]we show the ratio {Q\qq\Q) / {0\qq\0) eq. ()5.6p for pure neutron (left 
panel) and for symmetric nuclear matter (right panel). From the figure we observe that the leading order 
correction Si, eq. (|3.2p . is the dominant contribution. This gives rise to the linear density approximation 
represented by the dashed and dotted lines in the left panel and by the dashed line in the right one. The 
NLO corrections rising from Hg, eq. ()5.5p . for the case of neutron matter correspond to the difference 
between the solid and dashed lines and between the dot-dashed and dotted lines in the left panel of 
fig. [HI The former couple of lines refer to (0|nn|rj) and the latter to {Q\dd\Q). In this case the NLO 
corrections are small and amount to 4.5% of the LO ones at p = 0.3 fm^^. The corrections increase with 
density, as expected, since higher three-momenta are available for larger Fermi momentum. For the case 
of symmetric nuclear matter the NLO contributions are larger, though still mild. They correspond to the 
difference between the solid and dashed lines in the right panel of fig. [8] and they are 10% and 20% of the 
leading correction for p = 0.3 and 0.5 fm~^, respectively. The NLO corrections, Eq, tend to speed up the 
tendency towards a vanishing quark condensate in the nuclear medium (a signal of a possible restoration 
of chiral symmetry in nuclear matter). It is worth pointing out that the dependence on the subtraction 
constant go of the in-medium quark condensate to NLO is just at the level of a few per cent. This is 
shown in fig. [8] by the dot-dashed line where go = —0.5 m^ is used. For the reasons discussed, the quark 
condensate is significantly less dependent on go than E/A. Were cr = 64 MeV used |48j . the dominant 
linear density contribution will lead to a faster decrease of the in-medium quark condensate. E.g. at 
this level of approximation, the dashed line in the right panel of fig. [8] would cross the zero at already 
p « 2po. To this result one should add the difference between the solid and dashed lines corresponding to 
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the contributions from the in- medium nucleon-nucleon contributions, Hg, which are independent of the 
value of a taken. 

The one-pion exchange plus the contact nucleon-nucleon interaction terms from Cj^j^, eq. (j2.1ip . are 
fully iterated, as represented by the ellipsis in the diagrams 4-6 of fig. [TJ In ref. [30] one-pion exchange 
is iterated only once, and no nucleon-nucleon contact interactions at the same chiral order are included. 
However, the contribution of ref. [30] from once-iterated one-pion exchange is notoriously larger than ours. 
Keeping only this extra contribution, together with the linear density one, ref. [30] finds that the quark 
condensate vanishes already at p ~ 0.24 fm~^, while in our case this happens for the larger p = 0.41 fm~^ 
(this values coincides very closely to that for which E/A becomes positive in the right panel of fig. [7]) • This 
indicates that the further iteration of the one-pion exchange, together with the inclusion of the associated 
nucleon-nucleon local terms of /^tvjv ^^'^^ ^^ important impact. Other mechanisms were included in 



ref. [30] . They indicate a tendency towards a stabilization of the in-medium quark condensate in nuclear 
matter for densities above pQ. However, let us mention that no chiral power counting is followed by the 
authors in ref. [30] but an expansion in the number of loops. The latter does not coincide with a chiral 
power counting expansion because of the infrared enhancement associated with the nucleon propagators, 
as discussed above and also referred to in [51]. In addition, the local nucleon-nucleon interactions are not 
treated consistently within a chiral power counting either. Then, we consider as an interesting future task 
to include the higher orders needed within our power counting and confirm, if possible, the far reaching 
results of ref. [30] . Other papers also find that the linear decrease with density of the in- medium quark 
condensate is softened for p > po [271[32l[52] . However, important differences persist at the quantitative 
level. E.g. ref. [32] finds for densities above po a rnuch milder positive correction to the linear density 
decrease of the quark condensate than refs. |27 y 30 t [52]. Note that we obtain a realistic equation of state 
for symmetric nuclear matter and the quark condensate keeps its trend to vanish for higher densities, as 
also observed in ref. [32]. We agree on the observation performed in refs. [30^ 132] that the short-range 
nucleon-nucleon interactions effects are suppressed for the evaluation of the in-medium quark condensate 
which is dominated by long-ranged pion physics. This is a consequence of eq. (j4.13p . At the order we are 
working, if only the local nucleon-nucleon interactions from Cj^j^ were kept the derivative of Tjj with 
respect to m^ would be zero and Hg — t- 0. In the same way, the limit p^ — )• oo is independent of the 
pion mass which suppresses the influence of the short-range one-pion exchange part. It follows from the 
left panel in fig. [8] that for the case of pure neutron matter the linear tendency of the quark condensate 
is only weakly modified by the higher order corrections. Similar results are obtained in ref. [31] as well. 
Let us mention in passing that the approach of the Munich group [311153^ 154] has difficulties in providing 
a good reproduction of the equation of state for neutron matter while in our approach [2] it emerges in 
a quite straightforward way at NLO. 

6 Pion decay constant and the Gell-Mann— Oakes— Renner relation 

Let us now consider the axial- vector current A^ = g(x)7^75(rY2)g(j;), with q{x) a two-dimensional vector 
corresponding to the light quarks fields and r* the Pauli matrices. For notation we refer to ref. [21] . 

Due to the presence of the nuclear medium one should distinguish between the couplings of the pion 
to the spatial and temporal components of the axial-vector current, denoted by fs and ft, respectively. 
The fact that ft ^ fs in nuclear matter and some consequences thereof have already been discussed 
in [T7ti2n[55H59] . The calculations in |21j indicate a linear decrease of ft with density, ft = fn{^ — 
(0.26 ± 0.04)p/po)) where f-,^ = 92.4 MeV is the weak pion decay constant in vacuum and po is the 
nuclear matter saturation density. This result indicates that it makes sense to use chiral Lagrangians 
in the nuclear medium up to central nuclear densities. However, the result so far has been obtained by 
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Figure 9: Contributions to the in-medium pion decay up to NLO or 0{p^). The axial vector current is indicated 
by the wavy line and pions by the dashed ones. A wiggly line corresponds to the nucleon-nucleon interaction kernel 
(given in fig. [3]) whose iteration is denoted by the ellipsis. The diagram labelled tt-WFR indicates the contribution 
from the wave function renormalization of the pion. Crossed diagrams are not shown. 



considering pion-nucleon dynamics only, without including nucleon-nucleon interactions. 

We depict the contributions to the process of the in-medium pion decay (r2|yl;^|7r*) in fig. [H Diagrams 
1-3 were already considered in ref. [21]. Regarding diagram 4 one has to take into account the same 
comments as previously given in section [3] concerning diagrams 2 and 3. The contribution with only the 
free part of the nucleon propagators vanishes and that with only the density dependent part is taken 
into account by the diagrams 5 and 6. The remaining contributions, with one density dependent part 
for one nucleon propagator and a free one for the other, is suppressed by one order. The reason is 
because it implies the derivative of E^ with respect to energy. This already occurred for the diagram 2 
of fig. [1] in section [3l see also refs. [Il[2] for the case of the pion self-energy. For the diagrams 5 and 6 
we find a mutual cancellation. In the same way as we already found that for the case of the isovector 
contributions to the pion self-energy of diagrams (b) and (d) in fig. 4 of [Ij. Because of this there are 
no new contributions from the pion wave function renormalization (diagram indicated by vr-WFR in 
fig. [9]) beyond those already considered in ref. (21j . The driving mechanism for such cancellation has been 
explained here in section 14.11 for the case of the chiral quark condensate involving the mutual cancellation 
of S4 and H5. Due to suppression and cancellation of new contributions, we find that there are no 
additional contributions to those already given in [21] eq. (4.24) for the in-medium pion decay. 
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Together with the relation for the chiral quark condensate, eq. (j5.6p . we can write down the Gell- 
Mann-Oakes-Renner (GMOR) relation for symmetric nuclear matter [21] 



miff = -rh{n\uu + dd\n) + 60 , (6.1) 

where m-j^ is the in- medium pion mass and 60 corresponds to corrections that start at O(p^) in vacuum 
xPT, which implies that the GMOR relation is only exact at lowest order. The stability of the GMOR 
relation under the in- medium corrections as well as the fact that it is the temporal coupling ft, and not 
the spatial one fs, the one involved in the GMOR relation has previously been reported in |17H21| within 
the mean field approximation, and in |59j in the framework of QCD sum rules. Expanding the in-medium 
contributions to the pion mass, the pion decay constant and the chiral quark condensate explicitly we 
find 

mlf^ (1 + '^Sl + '^it + • • •) (^ + ^4? + 24? + •••)= -^^0|^^ + ^"^|0) (1 + 4'^ + 4'^ + •••)+ '^o , (6.2) 

where 62 , ^t and (5i denote the in-medium corrections of the pion mass squared, pion decay constant 
and chiral quark condensate, respectively, and the superscript indicates the corresponding chiral order. 
The relative 0{p'^) corrections on the left-hand-side of eq. (j6.2p . 5 { and 5r , only appear for the charged 
pions and correspond to isospin breaking due to different proton and neutron densities in asymmetric 
nuclear matter [21J . These in- medium corrections to the pion mass and the pion decay constant vanish 
for symmetric nuclear matter. They also cancel if we take the average of the pion masses or accordingly 
the pion decay constant. For these reasons, we can neglect the relative corrections at 0{p'^) for the 
consideration of the GMOR relation. As previously discussed, up to next-to-leading order the leading 

(3) 

nucleon-nucleon contributions to the pion decay, and therefore to the decay constant, cancel in 5\ . 
The same observation was already made for the pion self-energy in [l] and, therefore, contributions due 
to nucleon-nucleon interactions are absent in 6 2 as well. Regarding the right-hand-side of eq. (16. 2p . 

6^ has been calculated in section [3] and corresponds to Hi, eq. (|3.2p . In contrast to m^ and ft there 
is a non-vanishing NLO in-medium correction to the quark condensate due to the in-medium nucleon- 
nucleon interactions and given by Hg, eq. (I4.14p . Thus, 6^ 7^ 0. However, this is 0{p^), one order 
above the contributions discussed for the left-hand-side of eq. (j6.2p . Whence, it would be needed a full 
N^LO calculation for ft and ml in order to ascertain the stability of the in-medium corrections to the 
GMOR relation up to 0{p^), which is beyond the scope of the present work. On the other hand, as 
shown in [21j , the pion-nucleon dynamics that gives rise to 5 2 j ^f and 6~ for symmetric nuclear 



matter does not violate the in-medium GMOR relation, eq. (j6.2p . Then, we conclude that the in- medium 
corrections, including nucleon-nucleon interactions, do not spoil the validity of the Gell-Mann-Oakes- 
Renner (GMOR) relation up-to-and-including 0{p^) (or NLO) in our in-medium power counting scheme. 

7 Conclusions and outlook 

Employing the in-medium power counting of ^ and the methods of Unitary xPT for taking into ac- 
count the resummation of non-perturbative effects [2] we have calculated the chiral quark condensate 
up-to-and-including next-to-leading order, 0{p^) in nuclear matter. We have found an interesting partial 
cancellation between the diagrams involving the nucleon-nucleon interactions. In this way, those contri- 
butions that arise due to the leading quark mass dependence of the nucleon mass mutually cancel. This 
corresponds to the diagrams 4 and 5 in fig. [TJ As a result, only the diagrams 6 in the figure, that stem 
from the quark mass dependence of the pion mass, survive. This is the reason why previous calculations 
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have found that short range nucleon-nucleon interactions are suppressed for the calculation of the in- 
medium quark condensates, which is then dominated by the long-range pion contributions. Note that we 
did not exploit the Hellmann-Feynman theorem, but obtained our results explicitly from the generating 
functional, eq. (12. 2p . We conclude that the corrections are small for the quark condensate in the case of 
pure neutron matter. These corrections are more significant for the the case of symmetric nuclear matter. 
It is also worth pointing out that the full iteration of the nucleon-nucleon interactions reduces the force 
of such extra damping of the quark condensate as compared with other references. The dependence of 
our results on the subtraction constant qq is at the level of a few percent, much smaller than for the case 
of the binding energy [2j. On top of these nuclear effects a reduction of the present uncertainty in the 
pion-nucleon sigma-term would be most welcome. 

We have also addressed the calculations of the pion decay constant and tested the stability of the 
GMOR relation against in-medium NLO in our chiral counting, that is, up to 0{p^). The nucleon- 
nucleon contributions for the in- medium corrections vanish at this order, not only for the calculation 
of the pion mass but also for the in-medium pion decay constant. This implies that the expressions in 
ref. [21] for the coupling of the pion to the axial- vector current are right up to NLO, despite having 
considered the inclusion of nucleon-nucleon corrections. Because of this result we also find that in- 
medium contributions do not spoil the GMOR relation up to 0{p^), as in ref. [21]. The worked out NLO 
contributions to the quark condensate give rise to a new contribution to the GMOR relation at 0{p^) due 
to the nucleon-nucleon interactions that would require a N^LO calculation for the pion mass and weak 
pion decay constant in symmetric nuclear matter in order to check the stability of the GMOR relation. 
This would be an interesting future task. 

Higher order calculations are of most interest as they would provide the important two-pion exchange 
and multi-nucleon forces. They would furthermore merge meson-baryon mechanisms with novel multi- 
nucleon contributions that can be worked out systematically within our EFT. The large impact of the 
A-isobar excitation to symmetric nuclear matter in |30j, which is then partially compensated by the 
three-body interactions terms proportional to the low-energy constant ci, freezes the dropping of the in- 
medium chiral quark condensate at ~ 2pQ. This interesting fact requires confirmation within our power 
counting so as to keep all the terms contributing at the same chiral order, while keeping the full iteration 
of lowest order local and one-pion exchange diagrams. As commented above, the latter has a significant 
impact in the contribution at NLO. In addition, one should keep in mind that the proper way to address 
the issue of chiral symmetry restoration in the nuclear medium is the calculation of the temporal pion 
decay constant in the nuclear medium [21], which should be pursued for higher orders. 
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